Level structure and spin-orbit effects in semiconductor nanorod dots 
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We investigate theoretically how the spin-orbit Dresselhaus and Rashba effects influence the elec- 
tronic structure of quasi-one-dimensional semiconductor quantum dots, similar to those that can 
be formed inside semiconductor nanorods. We calculate electronic energy levels, eigen-functions, 
and effective gr-factors for coupled, double dots made out of different materials, especially GaAs 
and InSb. We show that by choosing the form of the lateral confinement, the contributions of the 
Dresselhaus and Rashba terms can be tuned and suppressed, and we consider several possible cases 
of interest. We also study how, by varying the parameters of the double- well confinement in the 
longitudinal direction, the effective gr-factor can be controlled to a large extent. 
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I. INTRODUCTION 

In recent years, much of the research in semiconduc- 
tor physics has been shifting towards spintronics^^ the 
novel branch of electronics in which the information is 
carried, at least in part, by the spin of the electrons. 
The electron spin might be used in the future to build 
quantum computing devices combining logic and storage 
based on spin-dependent effects in semiconductors. In 
order to achieve this goal, much study has been devoted 
recently to magnetic and optical 3 properties of semicon- 
ductors quantum dots££ and quantum wells£ One of the 
most popular spin-based devices was proposed by Datta 
and Das. 7 Improvements to the original design have been 
proposed recently by Egues et alr&S The Datta-Das de- 
vice makes use of the Rashba spin-orbit couplingi2*ii 
in order to perform controlled rotations of the spins of 
electrons passing through the channel of a field-effect 
transistor (FET), thus creating a spin-FET. The Rashba 
term is the manifestation of the spin-orbit interaction in 
quasi-one-dimensional (quasi-lD) semiconductor nanos- 
tructures lacking structural inversion symmetry. Addi- 
tionally, the lack of bulk inversion symmetry enables an- 
other spin-orbit term in the electronic Hamiltonian, the 
Dresselhaus termjiS which is also taken into account in 
the spin-FET design introduced in Ref. 

The influence of the Rashba and Dresselhaus Hamilto- 
nians in quantum dots (QD) has recently been treated 
in a number of theoretical works. The most-often 
studied geometry is that of quasi-two-dimensional dots 
with parabolic confinement in the plane lif^^i^ On the 
other hand, there is a growing interest and experimental 
progress in another type of quantum dots defined inside 
quasi-lD structures called nanorods or nanowhiskers^ 
In these structures, additional confinement in the lon- 
gitudinal direction can be introduced with great preci- 
sion, thereby allowing the formation of quasi-lD het- 
erostructures, such as multiple quantum dotsi 7 ^ and dot 



superlatticesi^S Nanorods can be grown out of numer- 
ous semiconductor materials. Their lateral widths can 
be controlled by selecting the size of the gold nanoparti- 
cles which are used to catalyze their growth and can be 
made as small as 3 nmiSi Recently, the transport prop- 
erties of these nanorod dots have been measured and the 
gated control of the number of electrons in them has been 
demonstrated^ 

Motivated by this experimental progress, we study the- 
oretically the electronic structure of quasi-lD coupled 
double dots including spin-orbit effects. This type of dot 
systems has also attracted interest in the field of quantum 
control of orbital wave functions due to their simplicity 
and tunability 21 ' 22 ' 23,24 As we will see here they are also 
well-suited for applications involving control of the spin 
degrees of freedom since they allow a great deal of control 
over the Rashba and Dresselhaus Hamiltonians. In this 
paper we study the influence of the Rashba and Dressel- 
haus spin-orbit Hamiltonians on the electronic structure 
of quasi-lD QDs, akin to those formed in semiconductor 
nanorods. Our emphasis on the spin-orbit interaction is 
obviously motivated by the current widespread interest 
in developing spintronic applications, which require a de- 
tailed understanding of the dynamics of the spin degree 
of freedom in semiconductor nanostructures. 

Let us denote by x and y the two transversal and by z 
the longitudinal direction of a quasi-lD nanorod, and let 
us call V z (z) the confining potential that defines a pair of 
coupled QDs along the nanorod. The laterally-confining 
potentials V x (x) and V y {y) are crucial in the determina- 
tion of the Rashba and Dresselhaus Hamiltonians and we 
consider different combinations of these potentials which 
can arise in our elongated geometry. We calculate the 
energy spectra and the wave functions by exact numeri- 
cal diagonalization of the total Hamiltonian and analyze 
how the energy levels and the effective ^-factor change 
as the Rashba and Dresselhaus couplings are modulated 
by varying the lateral confining potentials. Furthermore, 
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we study the effect of varying the size of one of the dots 
and the width of the central barrier between them. Since 
the strength of the spin-orbit interaction varies greatly 
among semiconductor compounds, we look at several ma- 
terials such as GaAs, InSb, GaSb, and InAs. Finally, we 
investigate the effective spin (S z ) as a function of the 
strength of the Rashba-like term for all the eigenhmc- 
tions of InSb with two different geometries. 

The quantization along different directions results in 
peculiar spin-momentum dependence. This in turn re- 
sults in SO effects that depend strongly on the symme- 
tries of the lateral confinement potentials. As such, the 
observation of SO spin splittings, as we will see, is directly 
attributable to asymmetry of the confinement and pro- 
vides an interesting probe of built-in strain fields and/or 
unbalanced composition gradients. 

We organize the paper as follows. In Sec. [H] we intro- 
duce the effective one-dimensional Hamiltonian and list 
the simplified forms it takes depending on the choice of 
confinement potentials. In Sec. lIIII we present the results 
for the energy levels including either the Dresselhaus or 
the Rashba term. In Sec. IIVI we study the effective g- 
factor and the expectation value of the z-component of 
the spin as a function of the strength of the Rashba term 
for different semiconductors and eigenstates. In Sec. |V] 
we provide a discussion and conclusion. 



II. THE ONE-DIMENSIONAL HAMILTONIAN 

We start with the complete Hamiltonian for a three- 
dimensional semiconductor structure in the absence of 
magnetic field, 



H = ^— + V(r) +H D + H R , (1) 

2m* 

where m* is the conduction-band effective mass, p is 
the momentum, V(r) is the confinement potential, and 
Ho and Hr are the general Dresselhaus and Rashba 
Hamiltoniansi2& Here we follow the current practice of 
calling Rashba terms those spin-orbit contributions to 
the Hamiltonian that arise due to the structural inver- 
sion asymmetry of the nanostructure, as opposed to the 
Dresselhaus terms which come from the bulk inversion 
asymmetry of the III-V semiconductors. Integrating out 
the x and y variables, we obtain the following effective 
one-dimensional Hamiltonian: 



H id = ^ + V z (z) + H ldD + H ldR , (2) 
2m* 
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FIG. 1: Potential-energy profile and schematic drawing of 
two Al/InSb coupled nanorod quantum dots. For InSb-based 
systems we take a well height of 100 meV, and for Al/GaAs, 
220 meV. In this example, the QD width is 300 A and barrier 
width 30 A, with smoothly changing barriers over a width of 
a few angstroms. The drawings (a)-(c) illustrate the lateral 
confinement geometries described in the text. 

H ldD = ^{<T X (p X )((p 2 y )-p 2 z )-<7y(Py)(( P l)-pl) 

+ o z p z ((pl)-(pl))}, (3) 

rr 1R 5 f/dVyi 3V Z \ 

where CTj, i = x,y,z, are the Pauli matrices, H\dD is 
the one-dimensional Dresselhaus term, and Hijr is the 
Rashba-like term enabled by the inversion asymmetry of 
the laterally confining potentials V x and V y . jr and 7d 
are parameters that depend on the materials. The aver- 
ages (...) are taken over the lowest-energy wavefunctions 
of the laterally confining potentials as we assume small 
nanorod widths. In Table I we present the parameters 
used in our calculations for different semiconductor ma- 
terials. An example of the confining potential in the lon- 
gitudinal direction, V z (z), is shown in Fig.Q] along with 
a schematic drawing of the nanorod QDs. 

We now list four different possibilities for the confin- 
ing potentials V x {x) and V y (y), based on the degree of 
symmetry of the structure. The Dresselhaus and Rashba 
Hamiltonians simplify considerably due to the fact that, 
in the absence of a magnetic field, the eigenstates can be 
chosen real, and therefore, expectation values of the mo- 
mentum are zeroi2& The four cases are (see Fig. l(a)-(c) 
for schematic drawings of the potentials in the first three 
cases): 

(a) Circular: V x (x) , V y (y) have inversion symmetry 
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about the origin and are equal, V x (x) = V y (y): 

Hso = HidD + HidR. = 0, (5) 
yields no SO contributions. 

(b) V x (x) , V y (y) have no inversion symmetry but are 
equal, V x (x) = V v (y): 

Hso = Hud + HxdR=^ \^^JPz{o x - cry) , (6) 

so that only Rashba terms are present. 
I 



(c) Elliptical: V x (x) , V y (y) are inversion symmetric 
functions and different, V x (x) ^ V y (y): 

Hso = H ldD + H ldR = -£<JzPz ((pl) - (pD) , (7) 



results in only Dresselhaus terms. 

(d) V x (x) , V y (y) have no inversion symmetry and are 
different, V x (x) =/= V y (y): 



H SO = H ldD + H ldR = j^-VzPz ((pi) - (pl)) + J^Pz Ux - °y (jfe/j > ( 8 ) 
represents the most general case and both Rashba and Dres selhaus contributions are present. 



For the calculation of the effective ^-factor we intro- 
duce a weak magnetic field along the z-direction. The 
field is chosen small so that the x — y orbital wave func- 
tions are not perturbed significantly. Thus, we only add 
a Zeeman term to the Hamiltonian, Hz = ^g^Baz, 
where /is is the Bohr magneton, B is the magnetic field, 
and go is the electron g-factor as per Table I. To calculate 
the energy levels and eigenfunctions, we expand the total 
Hamiltonian on a basis set of 300 wave functions of the 
quantum box of size L, i.e. 4> n . s (z) — \J^sin(J^)\ (s), 

where x( s ) is the spin function, and diagonalize it nu- 
merically without further approximations. The size L of 
the box is such that the whole double-dot structure is 
enclosed in the box, including the barriers on the sides of 
the dots and as such is irrelevant in the final eigenstates. 
We should notice that the geometry of the dots that we 
study here includes widths of 2-5 nm, while the most 
common nanorod widths in experiments are of the order 
of tens of nm. However, as we mentioned above, there 
are no experimental limitations to reducing the nanorod 
width to values we consider here. Smaller widths allow us 
to explore the basic physics and control of electronic wave 
functions with only one relevant lateral energy sublevel. 
Moreover, notice that typical charge depletion induced 
by the free surfaces further reduces the effective width of 
the nanorods, making them more ID-like. A final com- 
ment is that the incorporation of additional transverse 
levels in the nanorod is straightforward, but results in 
systems of coupled differential equations. 



III. ENERGY LEVELS 

We present results for the energy levels in cases (b) and 
(c), i.e. with only Rashba and Dresselhaus terms present, 
respectively. The general case (d) does not present quali- 



TABLE I: Parameters for semiconductor*— 



Parameter 


GaAs 


GaSb 


InAs 


InSb 


* / 27 

m = m/m<^ 


0.067 


0.041 


0.0239 


0.013 


„, 26/ a2\ 

7fl— ( A ) 


5.33 


33 


110 


500 




24 


187 


130 


220 


go 


-0.44 


-7.8 


-15 


-51 



tatively different features from (b) or (c) and therefore we 
concentrate here on the simpler cases. For case (b) we fix 
the strength of the Rashba term by giving the structural 
electric field (^jjr). For case (c), we use as confining po- 
tentials in the lateral directions two harmonic-oscillator 
potentials with different frequencies: V q (q) — ^m*uj q q 2 , 
q = x, y. These potentials have associated characteristic 
lengths t q = y / h/m*u) q . 

In Fig. we plot the two lowest energy levels for the 
InSb QDs taking (§£) = 0.5 me V/ A for case (b), and 
i x = 50 A, l y = 20 A for case (c). The indices on the 
horizontal axis denote the inclusion of different terms in 
the Hamiltonian. The figure shows how the energy lev- 
els of Hq (indices 1 and 4) are changed by the inclusion 
of a Rashba contribution H\ d R (case (b), index 2), and 
of a Dresselhaus contribution Hi d p (case (c), index 5), 
without magnetic field. With a weak magnetic field we 
have total Hamiltonians Ho + Hi dR + Hz (index 3) and 
Hq + H\ d n + Hz (index 6). We have carried out analo- 
gous calculations for the semiconductors quoted in Table 
I and the results were qualitatively similar to the ones 
shown here. The main general conclusion is that the ef- 
fect of H ld R is always stronger than that of H\ d u for the 
chosen parameters, which are representative of possible 
experimental situations. We note that the Rashba and 
Dresselhaus terms do not remove the spin degeneracy (as 
expected from the Kramers degeneracy in the absence of 
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FIG. 2: Ground-state and first-excited-state energy levels 
of the InSb nanorod QDs shown in Fig. 1. We compare the 

p 2 

eigenenergies of (1,4) Ho = + V z (z) to those of (2) Ho + 
H ldR , (5) Ho + H ldD , (3) Ho + H ldR + H z , and (6) H + 
H ldD +H z . B = 0.2T. 



magnetic field) but that they simply shift the levels down- 
wards, the strength of the shifts being controlled by the 
parameters (fg-) for Rashba and £ x and £ y for Dressel- 
haus. For the parameters chosen here the Rashba shift 
is of the order of O.lmeV for InSb and 0.1/iV for GaAs 
while the Dresselhaus shift is of the order of O.OlmeV for 
InSb and O.Ol^V for GaAs. 

As can be seen in Fig. [21 the energy shift produced by 
the H\dR varies quadratically with the structural elec- 
tric field (^r)- In Fig. 0] we show how the energy levels 
vary in case (c) as a function of £ x for the two lowest- 
energy states for fixed £ y — 50 A. The functional de- 
pendence here is also parabolic. This suggests that the 
spin-orbit corrections to the energy levels could be cal- 
culated fairly accurately with second-order perturbation 
theory. We performed the second-order perturbative cal- 
culation in the case with Rashba Hamiltonian, with a 
small magnetic field applied (0.1 T) in order to work 
with non-degenerate perturbation theory. A compari- 
son between the exact and second-order energies shows, 
for example, a difference of 17% for (jf-) = 1.5meV/A, 
and increasing differences for larger Rashba fields, as ex- 
pected. These results agree qualitatively with those of 
Ref. 01 f° r quasi-2D circular dots, where differences of 
up to 30% between the results of exact calculations and 
of second-order perturbation theory have been found. 



IV. EFFECTIVE ^-FACTOR 

The small magnetic field B = O.lTz breaks the spin 
degeneracy of the ground state and allows the calculation 
of the effective ^-factor (<?*) as a function of (^r) (case 
(b)) for GaAs, InSb, InAs and GaSb. In the figures we 
report normalized ^-factors: 
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FIG. 3: Contribution of the Rashba term to the energy levels 
of InSb (a) and GaAs (b) QDs as a function of (f£). GS: 
Ground state, 1 and 2: first and second excited states, respec- 
tively. Notice effect is much smaller in GaAs (energy given in 
/^eV), as anticipated. 
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FIG. 4: Contribution of the Dresselhaus term to the energy 
levels of InSb as a function of l x for the ground state (GS) 
and the first excited state (1) for £ y = 50 A. Level splitting 
in GaAs is barely visible on the same scale as in InSb. 



where E\ and E 2 are the Zeeman-split ground-state lev- 
els. Figure shows the results for case (b) (i.e. with 
only Rashba contributions) as a function of The 
decreasing trend of g* is qualitatively similar for all the 
materials but the magnitude of this Rashba effect varies 
greatly among them. The decrease of the g* is strongest 
for InSb and weakest for GaAs. 

We now examine what happens to g* when one modi- 
fies the features of the longitudinal potential V z (z), such 
as the barrier width w and the size of the QDs (so far 
we have taken Lqdi = Lqjj 2 = 300 A). In Fig. Efa) 
we show g* for w — 30, 130, and 330 A as a function 
of (tjj). We increase the barrier width but reducing at 
the same time the sizes of the two QDs so that the to- 
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FIG. 5: Effect of the Rashba Hamiltonian on the effective g- 
factor. g* /go for the ground state for different semiconductors 
as a function of (^r). 
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FIG. 6: Normalized effective g-factor for the ground state of 
InSb structures with Rashba Hamiltonian. (a) For different 
barrier widths w = 30, 130, 330A. (b) For different sizes of the 
QDs. Asymmetric case: Lqdi = 100 A and Lqd2 = 500 A; 
symmetric case Lqdi = Lqd2 = 300 A. 



FIG. 7: Mean value of S z and effective ^-factor for InSb sys- 
tems with symmetric V z (z) (two equal dots with Lqdi = 
Lqd2 = 300 A), (a) (S z ) as a function of (■fj) for the 
four lowest-energy doublets (pairs of Zeeman-split states), (b) 
g* /go for the same states. 



B = 0.1T is included. As expected, {S z ) = ±i in the ab- 
sence of ( ) . Figure 7 shows the results for a symmetric 
structure with Lqdi = Lqd2 — 300 A and Fig. 8 for an 
asymmetric one with Lqdi = 100 A and Lqd2 = 500 A. 
The symmetric case shows a crossing in (S z ) (Fig. 7(a)) 
while the asymmetric one does not (Fig. 8(a)). Using 
this information we recalculate the effective g-factor for 
the first four pairs of eigenstates for the symmetric (Fig. 
7(b)) and asymmetric (Fig. 8(b)) structures. The effec- 
tive g-factor, given here by the difference in (S z ) values 
for every Zeeman pair, vanishes at the crossing of {S z ). 
This vanishing of g* is a potentially useful effect in spin- 
tronics applications, as it can be achieved as a function 
of the potentially adjustable Rashba parameter (^jr)- It 
is interesting to note how different spatial asymmetry, 
introduced by the confinement potential along z (i.e. dif- 
ferent size dots), has strong effect on g* , and results in a 
finite value even at large Rashba fields. 



V. CONCLUSIONS 



tal size of the structure remains constant at 630 A. We 
note that increasing w leads gradually to having two un- 
coupled QDs and to a stronger variation of g* . In Fig. 
|Hfb) we set w = 30 A and change the QDs' sizes. We 
take Lqdi — 100 A and Lqd2 = 500 A in one case, and 
Lqdi — Lqd2 — 300 A in the other. We observe here 
that the symmetric potential produces a stronger varia- 
tion of g* than the asymmetric one. 

We look at these symmetric and asymmetric structures 
in more detail, and calculate the expectation value (S z ) 
as a function of (^r) for InSb dots and for the four lowest 
pairs of states (Zeeman doublets). Again a magnetic field 



We have studied how the spin-orbit Rashba and Dres- 
selhaus terms modify the electronic structure of nanorod 
quasi-one-dimensional double quantum dots. We have 
solved the problem by numerical diagonalization of the 
total Hamiltonian for varying confining potentials, in the 
lateral as well as in the longitudinal directions. The main 
conclusions of our work are the following: 
(1) For our system, the Rashba and Dressclhaus Hamilto- 
nians shift downwards the energy levels but do not break 
the spin degeneracy of the electronic levels in the ab- 
sence of an external magnetic field (as prescribed by the 
Kramers degeneracy.) 
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FIG. 8: Same as Fig. 7 for asymmetric V z (z) with Lqdi = 
100 A and L QD2 = 500 A 

(2) The Rashba effects are in general stronger than the 
Dresselhaus effects, but the latter are not negligible in 
general either. 

(3) Changing the strength of the spin-orbit terms, which 



is done by changing the lateral confinement length l x or 
£ y in the case of Dresselhaus or the structural electric 
field in the case of Rashba, results in energy levels 

that vary nearly quadratically with the control parame- 
ter. This indicates that the SO corrections to the energy 
levels are close to the second-order corrections in pertur- 
bation theory. We verified this result by comparing the 
exact and the perturbatively calculated energies. 
(4) By changing the strength of the Rashba term, the 
size of the central barrier, and the size and symmetry of 
the two QDs, it is possible to control the value of the ef- 
fective g-factor, which determines the Zeeman splitting. 
In particular, it is possible to make the effective (/-factor 
equal to zero. 



Acknowledgments 

We acknowledge support from the CMSS Program 
at Ohio University, Proycctos UBACyT 2001-2003 and 
2004-2007, Fundacion Antorchas, ANPCyT grant 03- 
11609, and NSF-CONICET through a US- Argentina- 
Brazil collaboration grant NSF 0336431. P.I.T. is a re- 
searcher of CONICET. 



1 G. A. Prinz, Science 282, 1660 (1998). 

2 S. A. Wolf et al., Science 294, 1488 (2001). Semiconductor 
Spintronics and Quantum Computation D. D. Awschalom, 
D. Loss, N. Samarth (eds) (Springer, Berlin, 2002). 

3 D. Gammon and D. G. Steel, Physics Today, October 2002, 
pp. 36-41. 

4 M. B. Tavernier et al, Phys. Rev. B 68, 205305 (2003); 
90, 146801 (2003). 

5 D. Loss and D. P. DiVincenzo, Phy. Rev. A 57, 120 (1998). 

6 K. J. Goldammer et al, J. Vac. Sci. Technol. B 16(3), 
May/June (1998). 

7 S. Datta and M. Das, Appl. Phys. Lett. 56, 665 (1990). 

8 J. C. Egues, G. Burkard, and D. Loss, Appl. Phys. Lett. 
82, 2658 (2003). 

9 J. Schliemann, J. C. Egues, and D. Loss, Phys. Rev. Lett. 
90, 146801 (2003). 

10 E. I. Rashba, Sov. Phys. Solid State 2, 1109 (1960). 

11 Y. A. Bychkov and E. I. Rashba, JETP Lett. 39, 78 (1984); 
J. Phys. C 17, 6039 (1984). 

12 G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

13 M. Governale, Phys. Rev. Lett. 89, 206802 (2002). 

14 E. Tsitsishvili, G. S. Lozano, A. O. Gogolin, Phys. Rev. B 
70, 115316 (2004). 

15 C. F. Destefani, S. E. Ulloa, and G. E. Marques, Phys. 
Rev. B 69, 125302 (2004). 



16 C. M. Lieber, Nano Lett. 2, 81 (2002). 

17 B. J. Ohlsson et al. Appl. Phys. Lett. 79, 3335 (2001). M. 
T. Bjork et al. Nano Lett. 2, 87 (2002). 

18 M. T. Bjork et al. Nano Lett. 4 1621 (2004). 

19 Y. Wu, R. Fan, and P. Yang, Nano Lett. 2, 83 (2002). 

20 A. M. Morales and C. M. Lieber, Science 279, 208 (1998). 

21 P. I. Tamborenea and H. Metiu, Phys. Rev. Lett. 83, 3912 
(1999); P. I. Tamborenea and H. Metiu, Europhys. Lett. 
53, 776 (2001). 

22 P. Zhang and X.-G. Zhao, Phys. Lett. A 271 419 (2000). 

23 L. Meza-Montes, S. E. Ulloa, and D. Pfannkuche, Phys- 
ica E 1, 274 (1998); A. Fendrik, M. J. Sanchez, and P. I. 
Tamborenea, Phys. Rev. B 63, 115313 (2001). 

24 C. E. Creffield and G. Platero, Phys. Rev. B 69, 165312 
(2004); Phys. Rev. B 65, 113304 (2002). 

25 C. F. Destefani, S. E. Ulloa, and G. E. Marques, Phys. 
Rev. B, accepted for publication. 

26 Principles of Quantum Mechanics, R. Shankar (Second 
Edition, Plenum Press, NY, 1987). 

27 O. Voskoboynikov, C. P. Lee, and O. Tretyak, Phy. Rev. 
B 63, 165306 (2001). 

28 V. I. Perel', S. A. Tarasenko, I. N. Yassievich, S. D. 
Ganichev, V. V. Bel'kov, and W. Prettl, Phys. Rev. B 
67, 201304 (2003). 



